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Bats are a unique group of mammals well suited to be hosts for emerging viruses. 
With current rates of deforestation and urbanization, redistribution of bat habitats 
to urban and suburban areas may bring bats into closer contact with livestock and 


humans. Common flying fox, Pteropus medius (previously known as Pteropus gigan- 


teus), forms large communal roosts on treetops, often in close proximity to human 
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habitation in Sri Lanka. This report describes the detection of coronavirus RNA in 
P. medius bat guano collected in Peradeniya, Sri Lanka. These viruses had >97% 


nucleotide identity with coronaviruses detected in Cynopterus sphinx, Scotophilus 


heathii and S. kuhlii bats in Thailand. Pteropus medius is widespread in Asia and 
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1 | INTRODUCTION 


Bats contribute to maintain a balanced ecosystem through dissemi- 
nating plant seeds or eating some types of insects. Some bat species 
are at risk and endangered. The displacement and redistribution of 
bats from their natural habitats expose susceptible human and ani- 
mal populations to micro-organisms (Han et al., 2015). Bats are a 
unique group of mammals that can fly long distances for daily forag- 
ing and aggregate in large colonies. Daily movement patterns, com- 
pact roosting structure, hibernation, periodic migration and long 
lifespan may make bats suitable hosts for emerging viruses (Han 
et al., 2015). 

Coronaviruses (CoV) are a group of enveloped, single-stranded 
(positive-sense) RNA viruses that belong to the order Nidovirales, 
family Coronaviridae. They are divided into four major groups: alpha, 
beta, gamma and delta coronaviruses, only the alpha and beta coron- 
aviruses being known to have infected humans (de Groot et al., 
2013). The earliest known human coronaviruses (HCoV), HCoV-229E 
and HCoV-OC43 are causes of mild upper respiratory infections. 


The identification of the precursor of Severe Acute Respiratory 


appears to excrete group D coronaviruses, which are hitherto confined to bats; 


however, these findings may have public health implications in the future. 


coronavirus, flying fox, Sri Lanka 


Syndrome coronavirus (SARS-CoV) in Rhinolophid bats led to more 
intensive investigations of coronaviruses in humans and in bats (Han 
et al., 2015). NL63 and HKU1 are newly discovered coronaviruses 
endemic in humans (Woo et al., 2005; Pyrc, Berkhout, & Hoek, 
2007), and MERS-coronavirus is a newly discovered zoonotic coron- 
avirus (de Groot et al., 2013). Recent studies have shown that bat 
coronaviruses were the source of many human coronaviruses includ- 
ing SARS-CoV, HCoV-229E and NL63 (Han et al., 2015). 


2 | MATERIALS AND METHODS 


A total of 51 treetops inhabited by fruit bats (Pteropus medius) were 
numbered on the first site visit, and these treetops were clustered in 
the Western corner of the Royal Botanical Gardens, bordering the 
River Mahaweli. A total of 50 bat guano samples were collected 
from randomly chosen treetops based on the original numbering. Of 
the 51 treetops numbered, 10 were chosen in the following order: 
1st, 5th, 10th, 15th, 20th, 25th, 30th, 35th, 40th and 45th and sam- 
pled from February to March 2014. The ground areas under the 
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selected trees were covered with clean polythene sheets. The accu- 





mulations of droppings were observed and identified immediately 
prior to collection to confirm the origin of the guano. About 1-2 g 
of each bat guano sample was collected into a vial with viral trans- 
port medium (VTM). The guano samples were transported to the 
Virology Laboratory, Faculty of Medicine, University of Peradeniya 
for screening with RT-PCR work. The ethical clearance (EC No: 
2012/EC/52) for the study was obtained from the Ethical Review 
Committee of the Faculty of Medicine, University of Peradeniya, Sri 
Lanka. 

Viral RNA was extracted using QIAGEN QlAamp Viral RNA kit 
(QIAGEN, Hilden, Germany) as recommended by the manufacturer. 
The cDNA was generated from RNA using Superscript Ill reverse 
transcriptase (Invitrogen) with random hexamers. The cDNA was 
screened for the presence of coronavirus using a broad range hemi- 


nested RT-PCR targeting conserved regions of the RdRp gene for 
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detecting known and unknown coronaviruses as previously described 
(Hemida et al., 2017). Positive and negative controls were included 
for each run. RT-PCR products were analysed using agarose gel elec- 
trophoresis to identify a product size of 440 base pairs indicating 
the presence of known or unknown coronaviruses. DNA sequencing 
was carried out on four strongly positive RT-PCR products at School 


of Public Health, University of Hong Kong, Hong Kong. 


3 | RESULTS AND DISCUSSION 


The bats from which the guano samples taken were morphologically 
identified as Pteropus medius (previously known as P. giganteus) 
belonging to the family Pteropodedae. Eight (16%) of the 50 bat 
guano samples were positive for coronavirus; three of 10 collected 


on 1 February, four of 10 collected on 8 February, zero of 16 


BtCoV Kenya Eidolon helvum/BtKY90 (GU065434.1) 
BtCoV Kenya Eidolon helvum/BtKY 102 (GU065376.1) 
BtCoV Kenya Eidolon helvum/BtKY50 (GU065394. 1) 
BtCoV Nigeria Eidolon helvum/13RS452-71 (KU131215.1) 
BtCoV Kenya Eidolon helvum/KY24/2006 (HQ728482.1) 
BtCoV Kenya Eidolon helvum/BtKY92 (GU065436.1) 
BtCoV Ghana Eidolon helvum/GH229 (KP231510.1) 
BtCoV Rousettus/HKU9-1 (NC_009021.1) 
0.88 BtCoV Rousettus aegyptiacus/Kenya/KY 06/2006 (HQ728483.1) 
BtCoV Rousettus/HKU9-2 (EF065514.1) 
BtCoV Dobsonia moluccensis/IFB2012-17F (AB918718.1) D 
0.97 BtCoV Rousettus/HKU9-5-1 (HM211098.1) 
BtCoV Ptenochirus jagori/2265/Philippines/2010 (AB683970.1) 
BtCoV BRT55629/Hipposideros lekaguli/CK/Tha/05/2012 (KJO20622.1) 
BtCoV Rousettus leschenaulti\GCCDC1 356 (KU762338.1) 
BtCoV Eonycteris spelaea/Singapore/2014 (KX452689.1) 
0.93 BtCoV B55700-1/Cynopterus sphinx/CB/Tha/5/2012 (KJ868721.1) 


BtCoV B55700-2/Cynopterus sphinx/CB/Tha/5/2012 (KJ868722.1) 
BtCoV B55700-3/Cynopterus sphinx/CB/Tha/5/2012 (KJO20604.1) 


Betacoronaviruses 


BtCoV Pteropus medius/Peradeniya-HKU CN12 
BtCoV Pteropus medius/Peradeniya-HKU CN15 
1 BtCoV B55736/Cynopterus sphinx/CB/Tha/6/2012 (KJO20606.1) 
BtCoV B55762/Scotophilus heathii/CB/Tha/6/2012 (KJO20607.1) 
BtCoV Pteropus medius/Peradeniya-HKU CN14 
BtCoV Pteropus medius/Peradeniya-HKU CN10 


BtCoV Rhinolophus/A1018/2005 (DQ648795.1) 
1— BtCoV Rhinolophus/HKU3 (DQ022305.2) 
HCoV SARS (NC 004718.3) 


BCoV (NC003045.1) 


0.82 HCoV OC43 (AY391777.1) 


MHV (FJ647223.1) 
HCoV HKU1 (NC 006577.2) 


wo 


> 


BtCoV Hipposideros/GhanaBoo/348/2008 (FJ710043.1) 


BtCoV HKU5-1 Pipistrelle/LMHO03f (NC009020.1) 


0.87 HCoV EMC (NC019843.3) 


IKE: BtCoV HKU4 Tylonycteris/HKU4-4 (DQ249241.1) 
BtCoV HKU6 Myotis pilosus/HKU6-1(DQ249224.1) 


0.98 


O 


BtCoV Myotis daubentonii/Germany/D8_ 38/2007 (EU375874.1) 


PEDV (NC 003436.1) 


0.81 BtCoV Myotis occultus/RockMountain27 (EF544564.2) 
BtCoV HKU2-1 Rhinolophus sinicus (DQ249235.1) 


HCoV 229E (NC002645.1) 
0.90 HCoV NL63 (NC 005831.2) 


Alphacoronaviruses 


BtCoV Miniopterus/AFCD307/03/06 (EU420137.1) 
0.95 BtCoV 1A Miniopterus/AFCD62/08/05 (EU420138.1) 


1 O81 BtCoV HKU7 Miniopterus (DQ249226.1) 
BtCoV HKU8 Miniopterus (NC 010438.1) 


7 FIPV (NC002306.3) 
TGEV (AJ271965.2) 


0.2 


FIGURE 1 Phylogenetic analysis of coronavirus RdRp gene nucleotide sequences. The tree was constructed by maximum-likelihood method 
using PhyML with an alignment of 339 nt. And the branch support was estimated by aLRT method, and values larger than 0.8 were indicated 
in the tree. The tree was mid-point rooted at the branch between alpha and betacoronaviruses. Bat coronaviruses from Sri Lanka identified in 
this study were highlighted in red. The groups A-D within the betacoronavirus lineage are indicated. GenBank accession number of sequences 


is shown in brackets 
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collected on 10 March and one of 14 collected on 19 March. Gen- 
ome sequencing of the RT-PCR products was carried out on four 
selected bat guano; BtCoV Peradeniya-HKU CN10 collected on 1 
February and BtCoV Peradeniya-HKU CN12, CN14 and CN15 col- 
lected on 8 February. The virus sequences showed >97% nucleotide 
identity with lineage D beta coronaviruses previously identified in 
Thailand in a fruit bat species, Cynopterus sphinx and insectivorous 
bats, Scotophilus heathii and Scotophilus kuhlii (Wacharapluesadee 
et al., 2015) (Figure 1). These viruses had approximately 75% 
nucleotide identity to the HKU9 like bat coronaviruses, which have 
widespread geographic distribution. 

Fruit bats, generally known as common flying foxes of genus 
Pteropus and family Pteropodidae, are a routine sight, in both rural 
and semi-urbanized habitats, in Sri Lanka. They may roost in trees 
closer to human habitation. In 2002, there were 30 bat species doc- 
umented in Sri Lanka (Yapa, Randeniya, & Ratnasooriya, 2002). 
Recently, a novel bat lyssavirus designated as Gannoruwa bat lyssa- 
virus was discovered in Pteropus medius bats in the same geographic 
area of Peradeniya (Gunawardena et al., 2016). Over the last few 
years, identification of new coronaviruses in bats including the most 
recently reported lineage D betacoronavirus in cave nectar bats 
(Eonycteris spelaea) in Singapore (Mendenhall et al., 2017) led us to 
investigate coronaviruses in fruit bats in Sri Lanka. Detection of lin- 
eage D betacoronavirus in guano collected from Pteropus medius bats 
shows the presence of these viruses in bats, which roost in public 
places such as the Botanic gardens. Lineage D betacoronaviruses are 
strictly confined to bats up to now; however, they might evolve to 
become harmful to humans with public health implications. Pteropus 
medius are one of the largest species of bats, roost in large trees 
near waterways or rivers and fly distances of up to 150 Km daily for 
foraging (Newman, Field, Jong, & Epstein, 2011). They are wide- 
spread in South Asia and can be commonly found in peri-domestic 
settings. In some cultures, these species are also hunted and con- 
sumed as food. New bats are born during February—May, and the 
samples collected in this study were from this period of new bat 
birthing and may explain the high rates of virus detection observed. 

A 7-week period of sampling has not enabled us to study the 
seasonality of the virus activity in the current study. The short frag- 
ment of the virus sequenced limits us from assigning a definitive tax- 
onomic status to what appears to be a new coronavirus species and 
the public health implications of these findings need to be further 


explored. 
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